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Abstract 

This article is focused on the dynamics of a rotating electrically conducting fluid in a turbulent 
state. As inside the Earth's core or in various industrial processes, a flow is altered by the presence 
of both background rotation and a large scale magnetic field. In this context, we present a set 
of three-dimensional direct numerical simulations of incompressible decaying turbulence. We focus 
on parameters similar to the ones encountered in geophysical and astrophysical flows, so that the 
Rossby number is small, the interaction parameter is large, but the Elsasser number, defining the 
ratio between Coriolis and Lorentz forces, is about unity. These simulations allow to quantify the 
effect of rotation and thus inertial waves on the growth of magnetic fluctuations due to Alfven 
waves. Rotation prevents the occurence of equipartition between kinetic and magnetic energies, with 
a reduction of magnetic energy at decreasing Elsasser number A. It also causes a decrease of energy 
transfer mediated by cubic correlations. In terms of flow structure, a decrease of A corresponds to 
an increase in the misalignment of velocity and magnetic field. 



1 Introduction 



Turbulence in electrically conducting fluids plays a role in many industrial processes and in flows from 
geophysical to astrophysical interests. In addition to containing a wide distribution of spatial and tempo- 
ral scales, these flows are often strongly anisotropic due to the effect of solid-body rotation, of a buoyancy 
gradient or of an imposed magnetic field. Independently, these forces and their effects on the anisotropy 
of turbulence have received an increasing interest. 



The case of rotating turbulence of non-conducting fluid has been extensively studied ( Jacquin et al. 
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homogeneous rotating turbulence has led to many advances in atmospheric and geophysical flows. A quick 
description of the main results can be drawn as follows (for more details, see Sagaut and Cambon 2008): 
the energy cascade is inhibited by the Coriolis force, so that the dissipation rate is reduced (Bardina 



et al. 1985 Jacquin et al. 1990); a transition from 3D to quasi-2D state develops and coherent vortical 
structures elongated in the direction of the rotation axis break the cyclone-anticyclone symmetry. 

A large scale magnetic field exists in many astrophysical objects (such as star and planet interiors or 
solar wind) or industrial configurations. The temporal and spatial variations of this imposed field may 



be neglected so that, at least locally, one may assume that it is uniform and steady (Moffatt 19711 



The transformation of initially isotropic turbulence by a static magnetic field is also a well documented 
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. The key phenomena are: a strong reduction 



of nonlinear transfers of energy along the direction of the imposed magnetic field, eventually leading to 
a quasi-twodimensional state and depending on the value of the imposed magnetic field; and anisotropic 
Ohmic dissipation whose relative amplitude depends on the magnetic Reynolds number. 
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In realistic flows, more than one force is usually responsible for anisotropic phenomena. The case 
of rotating stratified turbulence has been studied because of its direct application to atmospheric and 
oceanic flows (Cambon 2001 Liechtenstein et al. 20051. However, it seems that the coupled effect of 
the Coriolis force and the Lorentz force on homogeneous turbulence has received too few interest, except 
for theoretical considerations (Lehnert 1954 1955 Moffatt 1970). 

In this paper, we present the results of numerical simulations of incompressible homogeneous mag- 
nctohydrodynamic (MHD) turbulence submitted to a uniform magnetic field and to solid-body rotation. 
In order to avoid introducing additional parameters to the already complex parameter space (Reynolds 
numbers, Rossby number, Elsasser number, interaction parameter, spectral cut-off, etc), we consider the 
freely decaying turbulent regime, that is without introducing forcing terms. We also consider the simplest 
configuration, in which the imposed magnetic field is aligned with the rotation axis so that the config- 
uration is axisymmetric. We also select dimensionless numbers that correspond to regimes comparable 
to the ones observed in many geophysical and astrophysical systems: the Rossby number is very low 
(i.e. the rotation has a strong effect on the dynamics), the interaction parameter is moderate (i.e. the 
Lorentz force also has an important effect, but the regime is however far from wave turbulence), and the 
Elsasser number, characterizing the relative strength of the Coriolis and Lorentz forces, is about unity. 
Note however that the Lehnert number (see below) is always small compare to unity so that we consider 
flows dominated by rotation and inertial waves. 

In the following, we present the equations and nondimcnsional parameters in section 2, then the 
numerical resolution method in section 3, used to obtain the linear dynamics results presented in section 
4. Nonlinear simulations are discussed in section 5 starting with the characterization or rotating MHD 
turbulence in physical space, followed by an extensive discussion of the dynamics in spectral space. 
Conclusions are drawn in section 6. 



2 Governing equations and dimensionless numbers 

We consider the homogeneous turbulent flow of an incompressible electrically conducting fluid. The 
fluid is characterized by its kinematic viscosity v and its magnetic diffusivity rj; with r/ = (er/io)^ 1 , 
a the electrical conductivity and /j,q the magnetic permeability. The Navier-Stokes equations for an 
incompressible fluid, including the Coriolis and Lorentz forces, arc written in the rotating frame as: 

du 

— -vV 2 u = ux (u> + 2fl)+j x (B Q + b)- VP (1) 
at 

V-u = (2) 

where u> = V x u is the vorticity, b is the fluctuating magnetic field, j = V X b is the normalized 
electrical current and P is the pressure modified by magnetic pressure and centrifugal terms. Within the 
magnetohydrodynamic theory, the fluctuating magnetic field is derived from the induction equation 

db 

— - rfJ 2 b = Vx(ux(B Q + b)) (3) 

with the additional constraint 

V-6 = . (4) 

The initial state is fully developed and isotropic turbulence without magnetic fluctuations. The initial 
root mean square velocity is uq and the initial integral length scale is lo. The Reynolds number Re and 
its magnetic counterpart Rm are defined by 

Re = and R^i = . (5) 

V 71 

We focus on developed turbulent velocity field at large Reynolds number. The magnetic Reynolds number 
is also assumed large, with a magnetic Prandtl number of about one for all our simulations. Cases ranging 
from moderate to low magnetic Reynolds number, of specific interest for industrial applications, will be 
considered in future studies. 

The flow is submitted to a uniform magnetic field Bq (all the magnetic quantities are given in Alfven- 
speed unit such that B = B/^/p fj, where po is the density of the fluid) in a frame rotating with a 
constant angular velocity f2. Both Bq and $7 are vertical. The present axisymmetric case is the simplest 
configuration since all statistical quantities depend only on their orientations with respect to the direction 
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Table 1: Set of parameters of the initial condition for anisotropic DNS. R\ is the Reynolds number based 
on the Taylor microscale. 

Resolution up lp v R\ fcmax^ 

256 s 078 062 O0025 72 2.35 



of symmetry (the vertical one here). In the Earth's core, the strong differential rotation between the inner 
core and the mantle causes the magnetic field to be mainly toroidal, thus locally perpendicular to the 
rotation axis. However, this perpendicular case is much more complex in terms of anisotropy and is not 
considered here. The ratio between the eddy turnover time lp/up and the magnetic damping time rj/B 2 
is the magnetic interaction number 

^plp 



N = 



rju 



(6) 



This parameter can also be seen as the ratio between the amplitude of the Lorentz force with respect to 
the advective term in equation 0. In this paper, the intensity of the imposed magnetic field remains 
constant. We focus on the dynamics at moderate interaction paramter (N = 12, see below). Note however 
that we still consider the so-called strong turbulence limit where Bp ~ up, so that the anisotropy resulting 
from the applied magnetic field is small. The weak turbulence limit Bq 3> up has been considered in the 
non-rotating case (Bigot et al. 2008 Alexakis et al. 2007), and one observes an inhibition of nonlinear 



kinetic energy transfers in the direction of the imposed magnetic field. The effect of the rotation on this 
weak turbulence state will be considered in future studies. 

The ratio between the eddy turnover time and the Coriolis parameter is the Rossby number 



Ro = 



up 

2ni Q 



(7) 



In order to assess the effect of the Coriolis force on strong MHD turbulence, we consider different values 
of the rotation rate but we focus on the dynamics dominated by rotation so that the Rossby number is 
always small compared to unity. The product between the above two dimensionless parameters defines 
the Elsasser number 

characterizing the respective influence of the Coriolis force and the Lorentz force. For non-rotating MHD 
turbulence, A — > oo and we consider in the following Elsasser numbers down to 0.5. Note that it is 
believed that the Elsasser number of the Earth's core is about unity. The relative importance of the 



rotation and of the magnetic field can also be quantified by the Lehnert number ( Lehnert 1955 ) , which 



is the ratio between the frequency of Alfven waves and the frequency of inertial waves: 



C 



Bp 
2M 



(9) 



In the present simulations, the Lehnert number is always small compared to one so that inertial waves 
are rapid compared to Alfven waves. To reach larger values of the Elsasser number, one has to increase 
the value of the imposed magnetic field, thus considering the weak turbulence limit. We also introduce a 
scale-dependent Lehnert number, based on a wave number k, as 



C = 



B k 
2fl 



(10) 



3 Numerical method 

A Fourier pseudo-spectral method is used to solve equations to Q. The velocity field and the 
fluctuating magnetic field are computed in a cubic box of side 2ir with periodic boundary conditions 
and a resolution of 256 3 Fourier modes. A spherical truncation of Fourier components is used to remove 
completely the aliasing and the time scheme is third-order Adams-Bashforth. The dissipative terms 
(proportional to v and 77 in ([T]) and ([3])) are implicitly solved through the change of variables u(k,t) <-» 
u(k, t) exp {—vk 2 i) and b(k, t) <-> b(k, t) exp (~i]k 2 t) where u(k, t) and b(k, t) are the Fourier transforms 
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Figure 1: (Color online) Craya-Herring frame (e^\ e^ 2 \ e^ 3 ') in Fourier space. An incompressible Fourier 
variable is perpendicular to the local wave vector k so that its contribution along e^ 3 ) is zero. Fourier 
modes in the shaded region contribute to E(k, 9) (see equation (25)). The polar modes (9 ~ 0) contribute 
to horizontal kinetic energy, whereas equatorial modes (9 ~ 7r/2) contribute to both vertical (along e^) 
and horizontal (along e^) energies. 



of u(x, t) and b(x, t) respectively. The nonlinear terms are computed in physical space, and then obtained 
in Fourier space by using fast Fourier transforms. 

The initial Eulcrian velocity field comes from an isotropic simulation in the purely hydrodynamic case 
with zero angular velocity, whose initial velocity field is a random superposition of Fourier modes with 



an energy distribution given by the narrow-band spectrum E(k,t = 0) ~ k exp(— 2(k/ki) 2 ) (Rogallo 



1981). During the initial stage only, a large-scale forcing is applied until a statistical steady state of 



classical isotropic turbulence is reached. This forcing is applied for wavenumbers such that 2 < k < 4. 
The resolution is such that the minimum value for k max l rj is 2.35 for all computations, where l v is the 



Kolmogorov length scale (Jimenez et al. 1993). At the beginning of the anisotropic simulations, the 



magnetic field Bq and the angular velocity are suddenly switched on, the forcing is turned off, and there 
are no magnetic fluctuations. We choose to consider the case without initial magnetic fluctuations to 
reduce the already wide parametric space and to focus on the growth of magnetic fluctuations. If the flow 
contains initial magnetic fluctuations, one has to define additional parameters, such as the initial ratio 
between magnetic and kinetic energies or velocity/magnetic field cross-correlations. Some characteristics 
of the initial condition used for anisotropic simulations are gathered in Table [T] 



4 Linear inviscid dynamics 



Let us first consider the inviscid linear dynamics of the problem. Some details about the following 



oscillating solutions can be found in Lehnert (19541 and Moffatt (1971). It is convenient to introduce the 
so-called Craya-Herring frame (see figure |1] and Sagaut and Cambon ( 2008 ) for the link to the toroidal- 
poloidal decomposition in physical space) defined by 



,(i) - 



A^L, e O) = fc(D and e <3) * 
\k x n\ k k 



(11) 



Given that both velocity and magnetic fields are divergence-free, the components along the wave vector 
are zero so that u^- 3 \k) = U 3 ^(k) = 0. Neglecting nonlinear terms and dissipative effects in equations 
and ([3]), one can write the following linear system in the Craya-Herring frame: 



( 



U>i 

V 



-Uli 




-iw„ 



-lW a 









-iw a 










6d) 
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w 



= 



(12) 
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Figure 2: Visualizations of magneto-inertial wave packets. The vertical component u z of the velocity is 
presented in a vertical plan containing the initial impulse. u™ ax is the maximum value of the vertical 
velocity. The dots indicate the vertical position of the initial perturbation. All visualizations are presented 
at tBo/lf fa 6 where If is the characteristic length scale of the initial perturbation, (a) Pure Alfven waves 
(B = 1 and fl = 0). (b) Magneto-inertial waves with B = 1 and = 3. (c) Magneto-inertial waves 
with Bq = 1 and SI = 10. (d) Pure inertial waves (Bq = and ft = 10). 



where u>i — 2Q cos 9 is the dispersion relation of inertial waves, iu a = B k cos 6 is the dispersion relation of 
Alfven waves and 9 is the polar angle between i~2 (or Bq) and the wave vector k. The dispersion relation 
of magneto-inertial waves follows from the eigenvalue problem, and write 

w = \ (±w< ± V^+*^) ■ ( 13 ) 

In the following, we assume that the Lehnert number is small so that w, 3> to ai i.e. the characteristic 
frequency of inertial waves is much larger than the one of Alfven waves. Note that this assumption is 



valid only for a certain range of scale since the Lehnert number defined by equation ( 10 1 depends on 
the wave number k. Using this assumption, and considering that the initial condition is u(k,t = 0) = 
(u^(k, 0), u^(k, 0)) for the velocity and b(k,t = 0) = (0,0) for the fluctuating magnetic field, the 
solutions are 

(uW(k,t)\ _ [axuit -aauA fu^(k,0)\ . . 

\u^{k,t)) \smwit cosuj.t J \u^(k,Q)J [ ' 

(U 1 \k,t)\_ 2w a. (i , /9) ( cosM/2) sin( Wl t/2)\ fu^(k,0)\ 

{bM(k,t)J ~V (i/j V- sin (^*/ 2 ) cosM/2)J U (2) (fc,0)J ' (15) 

The time dependence of the velocity Fourier mode is characteristic of pure inertial waves: the Fourier 
mode rotates in the plane (e^, e' 2 )) with an angular velocity fixed by the dispersion relation of inertial 
waves. Concerning the magnetic mode, we first note that its amplitude is modulated by uj a /uji, which is 
assumed to be very small here, so the magnetic energy should be very small compared to the kinetic energy 
when the rotation is dominant. Secondly, it is known that, without rotation, the magnetic mode and the 
velocity mode are colinear (see for example Moffatt 1967) so that pure Alfven waves are characterized 
by u(k,t) |j b(k,t) (independently of the wave vector considered, so that this preferential alignment still 
holds in physical space). With rotation, a phase shift appears between the magnetic and kinetic modes 
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Figure 3: (Color online) (a) Kinetic energy (continuous lines) and magnetic energy (dotted lines) for 
magneto-inertial waves. The time is normalized by lf/B where If is the characteristic lenght scale of 
the initial impulse and -Bo is the imposed magnetic field (in Alfven speed units). (b) Probability density 
function of the cross correlation (16) at tBo/lf ~ 10. 



in spectral space, so that there is no reason to observe velocity fluctuations parallel to magnetic ones in 
physical space. 

In the following, we propose a simple numerical experiment to assess qualitatively the previous linear 
statements. Using the numerical method described above, we perform direct simulations of magneto- 
inertial waves, with a resolution of 128 x 128 x 256 Fourier modes. This resolution is well above what 
is needed for such simulation but it ensures that no spurious effects of the periodic boundary conditions 
appear. From a non- magnetized fluid initially at rest, we introduce a small horizontal velocity perturba- 
tion (both Jl and Bq are vertical) localized in space and in time to approximate the impulse response of 
the system. The amplitude of the perturbation is very small so that nonlinear interactions are negligible 
throughout the simulation (nonlinear terms are nevertheless removed to ensure that the results are purely 
linear). We set Bq — 1 and we consider four rotation rates: S7 = (pure Alfven waves), 3, 10 and 20, and 
we also perform a purely hydrodynamic simulation with B = and f2 = 10. A volume rendering of the 
vertical component of the velocity is presented on figure [2] Without rotation, an upward and a downward 
Alfven wave packets propagate along the imposed vertical magnetic field at a constant velocity. The shape 
of each wave packets is directly linked to the initial perturbation, and does not vary in time (dissipative 
effects are neglected). When rotation is introduced, dispersion of the energy along non- vertical directions 
is observed, which is very similar to the one observed in purely inertial waves (see figure [2jd) ) . It is 
well-known that the frequency of inertial waves depends only the angle 9 between ft and k. In the case 



of a sinusoidal forcing (McEwan 1970), the ratio between the forcing frequency and the rotation rate 
sets the direction of propagation of energy. In our case, all the available frequencies are excited by the 
initial pulse so that one observes different directions of propagation. Figure |3][a) presents the evolution 
with time of both kinetic and magnetic energies. Without rotation, the equipartition between kinetic 



and magnetic energies is observed after some relaxation. As predicted by the linear solution (15), the 
magnetic energy is strongly damped in dominant rotation cases. Finally, the probability density function 
of the cross-correlation between u(x) and b(x), defined as 

2(ii(ai).6(aO) 

P\ x ) — x ,n/ \ (Id) 

v ; u 2 (x) + b 2 (x) v ' 

is plotted in figure |3jb). For pure Alfven waves, the colinearity between u and b is clear and the rotation 
tends to break the preferential alignment. For low value of the Lehnert number (i.e. dominant rotation), 
it seems that the perpendicular state becomes more probable. 

The previous linear observations are to be kept in mind for the following nonlinear study. However, 
in order to be clear, we insist on the fact that the following MHD direct simulations are in the strong 
turbulence regime, meaning that nonlinearities are dominant with respect to the imposed magnetic field. 
In that case, the preferential alignment between u and b will not be observed but the linear effect of the 
rotation could give some insights on the nonlinear results. 
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Table 2: Initial values of the parameters for DNS computations presented in section [5] 



Run 


v = r\ 
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B 


Rm 


Ro 


N 


A 
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0.0025 
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0.0025 
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00 
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0.0025 
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0.16 
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0.2 


200 
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0.0025 


16 


0.2 
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0.04 


12 


0.5 


0.01 



5 Nonlinear simulation results 

We now focus on the nonlinear turbulence regime of initially homogeneous isotropic developed turbulence 
submitted to both background rotation and a uniform steady magnetic field. In that case, the growth 
of magnetic fluctuations is due to the deformation of the uniform magnetic field lines by the velocity 
fluctuations, leading to the generation of Alfven waves. In all the runs, the magnetic fluctuations are 
initially zero. We perform different anisotropic simulations to understand the effect of rotation on the 
growth of magnetic fluctuations. The main parameters are gathered in Table [2| First, the value of the 
magnetic diffusivity is fixed so that the initial magnetic Reynolds number based on the integral scale 
is about 200, the magnetic Prandtl number being equal to 1. The value of the imposed magnetic field 
corresponds to an interaction parameter of about 12. We consider three different rotation rates (Runs 1 
to 3 in Table [2]), so that the Rossby number varies from 0.6 to 0.04. For reference, we also perform an 
isotropic hydrodynamic simulation (run ISO in Table [2]) and a magnetohydrodynamic simulation without 
rotation but with an imposed magnetic field (run MHD in Table [2]). At time t — 0, the magnetic field B a 
and the rotation Jl are applied. Note that the following results are qualitatively unchanged if the rotation 
is introduced before the magnetic field such that anisotropy due to the Coriolis force is already present 
when the magnetic field is applied. However, the effect of the Coriolis force on pre-existing magnetic 
fluctuations is a different matter and will be considered in future works. In the following, the time is 
non-dimensionalized by the initial eddy turnover time, t* = tuo/lo. 

5.1 Characterization in physical space 

First, we propose a quick characterization of the flow in physical space. Visualizations of the current 
density at time t* = 5 are gathered in figure [4] The non-rotating results are on the left whereas the 
rotating results are on the right. Without rotation, the current density is clearly concentrated in sheet- 




Figure 4: (Color online) Visualizations of the square current density j 2 at t* — 5. is about 25% of 
the maximum value in the numerical domain, (a) Non-rotating MHD case, (b) Rotating MHD case with 
Aw 0.5. 
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Figure 5: (Color online) (a) Normalized pdf of the correlation coefficient p(x) between u and b at t* = 5. 
(b) Normalized pdf of fo/ert,, where ab is the variance of b, at t* —5. The symbols correspond to horizontal 
components whereas solid lines correspond to vertical components. 



like structures randomly orientated. Even if the interaction parameter is relatively large, we do not 
observe well-defined anisotropy since the flow is still in a strong turbulence regime (i.e. B w u ). In the 
rotating case (A « 0.5 presented in figure |4|b)), the current sheets are clearly aligned with the rotation 
axis. Note in addition, that the current sheets are thinner in the rotating case than in the non-rotating 
case. Concerning the velocity field in the rotating case (not shown), the visualizations are very similar to 
what is observed in hydrodynamic rotating turbulence, since the imposed magnetic field is weak: vortices 
are elongated along the rotation axis. 

The correlation between the velocity and magnetic field fluctuations can be measured by the cross- 



correlation coefficient defined by equation (16). This quantity has already been used for the linear 



simulations of magneto-inertial waves presented in section [4j In the non-rotating case, more detailed 



results about this quantity in the context of MHD turbulence can be found in Bigot et al. (2008). First, 
figure [5ja) presents the probability density function of p(x) at the end of the computation (i.e. for 
t* =5). Independently of the rotation rate, the pdf are centered around zero (in fact, in the non-rotating 
case, the pdf is slightly shifted towards negative values, see Bigot et al. 2008 for more details). The 
extremal values p = — 1 and p — 1 correspond to downward and upward propagating Alfven waves. 
Note that, in contrast with linear prediction, we do not observe a dominance of Alfvenic fluctuations. 
However, in the strong turbulence regime we consider here, nonlinearities are dominant so that the 
preferential alignment of Alfven waves is broken. As the rotation increases, one observes that the pdf 
peaks around zero, indicating that the magnetic fluctuations are mainly oriented perpendicularly to the 
local velocity fluctuations. This is consistent with the previous linear analysis. The large scale fluctuations 
are dominated by inertial waves so that Alfvenic fluctuations (characterized by an equipartition between 
kinetic and magnetic energies, and by rs I) are damped. 

The pdf of the vertical (solid lines) and horizontal components (symbols) of the fluctuating magnetic 
field are gathered on figure U». One can sec that the horizontal component of the magnetic field is 
dominant and more intermittent in rotating cases than in non-rotating case, which is consistent with the 
visualizations, where very thin current sheets are aligned with the rotation axis. 

5.2 Dynamics 

Figure [6Ja) presents the time evolution of the magnetic (bottom of the figure) and kinetic energies (top 
of the figure) for different values of the Elsasser number. In addition, figure ^b) represents the ratio 
between the kinetic energy and the magnetic energy, also called Alfven ratio. In the wave turbulence 
regime (i.e. Bq 3> uq), this ratio is expected to be equal to one, since Alfven waves are dominant and 
characterized by an equipartition between kinetic and magnetic energies. However, non-rotating strong 
MHD turbulence submitted to a uniform magnetic field is characterized by an Alfven ratio slightly smaller 
than one, indicating the presence of non- Alfvenic fluctuations (see e.g. Bigot et al. 2008). In the non- 
rotating case (see the A — > oo curve on figure [6^b)), this quasi-equipartition state is indeed observed 
at the end of the simulation (i.e. t* > 4), with a slight dominance of magnetic energy. However, as 
the Elsasser number decreases, i.e. as the rotation rate increases whereas the imposed magnetic field 
remains constant, the magnetic fluctuations are strongly damped whereas the kinetic energy decrease is 
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slower than in non-rotating cases. Thus, the Alfven ratio increases with the rotation rate, indicating that 
Alfvenic fluctuations are dominated by another phenomenon, namely inertial waves. This is in accordance 
with the linear predictions of section[4j and the similarities between figures[3]ja) and[6^a) show that waves 
play an important part in the dynamics. The line with symbols on figure |6][a) corresponds to the isotropic 
computation (i.e. Bq = and £1 = 0) obtained from the same initial condition. Without rotation (i.e. 
A — > oo), the kinetic energy decays faster in the magnetized case than in the hydrodynamic case due 
to the additional ohmic dissipation. As the rotation rate increases, the kinetic energy decay rate is 



slower, which is a well-known result already observed for non-magnetized rotating turbulence (Jacquin 



et al. 1990). However, the damping of the magnetic energy due to background rotation is an original 
observation. 

We define the kinetic dissipation rate tx and the magnetic dissipation rate €m as 



e K = v (uj 2 ) 



(17) 
(18) 



ex and em are plotted in figures [7|a) and (b) respectively. The kinetic dissipation rate is always decreas- 
ing, since we compute MHD flows from developed turbulence already containing small-scale fluctuations, 
and is always smaller than its isotropic value (line with symbols in figQa)). It seems inconsistent to 
observe a smaller value of the kinetic dissipation between isotropic and MHD turbulence, since the latter 
contains additional dissipative effects. Note however that in the present MHD case without initial mag- 
netic fluctuations, some of the initial kinetic energy is transferred to the magnetic energy, thus reducing 
the amplitude of the kinetic motion and dissipation. Note in addition that the sum of both dissipation 
rates ex and cm (not shown) is greater than the isotropic kinetic dissipation rate. The initial effect of 



rotation is, as in non-magnetized cases (Cambon, 2001 Morinishi et al. 2001), to decrease the kinetic 



energy dissipation. However, at larger times (i.e. for t* > 2), ck is greater for rotating cases than for 
non-rotating ones. This result deserves some explanations. The well-known effect of the Coriolis force is 
to decrease kinetic energy transfer (see figure |8][a) below). This implies a reduction of ex at the begin- 
ning of the simulation. However, as time increases, the decay of is also reduced by rotation (since the 
kinetic energy is less dissipated), so that it can exceed its non- rotating value. 

Concerning the magnetic field, we observe an initial increase of the magnetic dissipation rate due 
to the generation of small-scale magnetic fluctuations. After reaching a maximum, €m decreases. As 
rotation is introduced, the small-scale development of the magnetic fluctuations is accelerated (since the 
maximum of ejv/ occurs at earlier times) but the maximum of the magnetic dissipation is substantially 
reduced. The reduction of cm with the rotation rate compensates the increase of ck at large times, so 
that the whole dissipation rate ex + cm always decreases with the Rossby number. 

Let us now discuss the transfers of both kinetic and magnetic energies. In the equation ([I]) governing 
u, in addition to the classical nonlinear advective term = u X u>, the nonlinear term which reflects 
the impact of the fluctuating magnetic field is = j x b = (V x b) x b. Similarly, the advective term 
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Figure 6: (Color online) (a) Kinetic energy K and magnetic energy M versus dimensionless time t* 
tu /l . (b) Alfven ratio K/M. 
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Figure 7: (Color online) (a) Kinetic energy dissipation rate ex versus time, (b) Magnetic energy dissi- 
pation rate e^- versus time. 



in the induction equation pi is g = V x (m x 6), which can be rewritten as 



d , , , , d „ . d 



11, ; 



n = i^r^ b 3 - U A0 = -^rihuj) + jr-bj (19) 



dxj dxj dxj 



with implicit summation over repeated indices. This equation can be recovered using Ricci formulae, but 
is well known with a complete analogy with the Helmholtz equation for fluctuating vorticity. 
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Figure 8: (Color online) (a) Kinetic energy transfer spectrum Tx(k). The dashed thick line corresponds 
to the initial transfer spectrum at t* = 0. Solid lines correspond to transfer spectra at t* = 2. (b) 
Total energy transfer spectrum Tx(k). Dashed lines correspond to transfer spectra at t* = 0.2 whereas 
solid lines correspond to t* —2. Given that the figures are plotted in semi-log scale, the transfers are 
multiplied by k to enhance the fact that their integrals are zero, (c) Spectra of b-q. The dashed line 
corresponds to the spectrum at t* = 0.2 whereas solid lines correspond to spectra at t* = 2. 
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The dot product of equation ([I]) by u yields the equation for the kinetic energy K = ^u-u, or 
dK 

— vu-Xl 2 u = M-s (tl) + 2w (u x n) + u-(j x B ) +m-s (6) - u-VP . (20) 

at ^_ 

(o) 



Terms which are obviously zero, such as u-s^ u \ are retained in equation (20) in order to anticipate the 
origin of "true" spectral transfer terms with zero integral. The dot product of equation ^ by b yields 
the equation for the magnetic energy M = \b-b, or 

^ -rjb-V 2 b = b-(V x (m x B Q ))+b-q . (21) 
at <• „ / 

(b) 

Important cubic terms in the kinetic energy equation are 

^'l^-t^-llf'^ < 22) 



Their counterpart in equation (21) for M are 



&-9=^M J -^(6V). (23) 
oxj 2 axi 



Energy conservation is found by taking ensemble averages of equations (20) and (21 1 and ignoring 



dissipation terms. Mean field terms play no role: the second term in the right hand side of equation 



( 20 ) is zero because the Coriolis force produces no energy whereas the terms (a) and (6) in equations 



(20) and (21) linked to Alfven velocity cancel each others. Only the cubic terms in equations (22) and 
(23) control the energy transfer. Terms which appear as the gradient of a cubic product are zero due to 
statistical homogeneity and their spectra really correspond to a true downscale transfer. In contrast, the 
term (dui/dxjbibj) is similar to the nonlinear vortex stretching term in the hydrodynamic case, replacing 
6 by w. It is nonzero and does not correspond to a transfer term. In counterpart, this term appears 



with opposite sign in the kinetic equation (20) and the magnetic equation (21), so that it is a pure 
'intercomponent' (from magnetic to kinetic) energy transfer, although not a true 'scale-to-scale' energy 
transfer, and has no contribution to the equation for total energy K + M. 

Figure |8ta) presents the classical spherically averaged transfer spectra Tx{k) derived from the ad- 



vective term u-s^ u ' in equation (20). We plot the value of the transfer multiplied by the wave number 
so that, in semi-log scale, one clearly sees that the positive area is equal to the negative one. As for 
non-magnetized rotating turbulence, the rotation tends to reduce nonlinear transfer from large scales 
to small scales, thus reducing the overall kinetic energy dissipation (at least initially, see figure (7^a)). 
Figure |8|b) presents the transfer spectra Tr(fc) derived from the sum m-s^ + b-q, which correspond to 
the spectral transfer in the equation for the total energy. Again, the effect of the rotation is clearly to 
reduce nonlinear transfers, thus reducing the total (i.e. kinetic plus magnetic) dissipation. Finally, the 
advection of the magnetic fluctuations by the velocity field is characterized by the term b-q. This term 
is not a transfer (its integral over k is non zero) but is however also reduced by the rotation, as shown 
on figure [8](c). 

Spherically averaged energy spectra are represented in figure [9j for four different Elsasser numbers 
from oo (non-rotating case) to 0.5. All the spectra are plotted at the same dimensionless time t* « 5. 
The magnetic energy spectra EM{k) correspond to the dotted blue lines whereas the kinetic energy spec- 
tra correspond to the solid red lines. In each graph, the initial kinetic energy spectrum is presented 
for reference as a thin solid line. We retrieve the global reduction of the magnetic energy while in- 
creasing the rotation rate, as already observed in figure [6j Without rotation (A — > oo, figure [9ja)), the 
kinetic/magnetic quasi-equipartition of energy is observed at all scales of the flow, with a slight excess of 
magnetic energy with respect to the kinetic one for inertial and dissipative scales. At moderate Elsasser 
number (A w 8, figure |9jb)), the rotation affects predominantly the large scales of the flow, i.e. the 
magnetic fluctuations are attenuated at small wave numbers, whereas the large wave numbers are still 
equipartitioned, with an Alfven ratio about unity. When rotation increases (figures [9^c) and (d)), the 
magnetic energy is damped over all scales of the flow, whereas the kinetic energy spectrum at small wave 
numbers recovers the initial isotropic level, which is consistent with the decrease of the kinetic direct 
cascade. This scale-dependent attenuation of the magnetic energy is directly linked to the competition 
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Figure 9: (Color online) Kinetic energy spectra Ex{k) and magnetic energy spectra EMik) at t* = 5. 
Vertical arrows corresponds to the cut-off wave number k c = 20,/ B . The fc -5 / 3 and k~ 3 slopes are 
shown for comparison with classical scalings in hydrodynamic isotropic and rotating turbulence (see for 
example Sagaut and Cambon 2008). (a) A -> oo. (b) A = 8. (c) A = 2. (d) A = 0.5. 



between Alfvenic fluctuations and inertial waves due to rotation, as shown by the previous linear analysis. 

It is important to note that the dispersion properties of inertial waves are very different from Alfven 
waves. At a given scale k, the maximum frequency of Alfven waves is B k whereas the maximum 
frequency of inertial waves is 2f2. By equating these two frequencies, one can therefore estimate a cut-off 
wave number k c = 20/ Bq. The value of k c is indicated by a vertical arrow in figures [9jb) and (c), which 
clearly separates two spectral domains in wavespace. In the spectral domain such that k < k c , inertial 
waves are the fastest phenomenon, so that the energy transfers from kinetic to magnetic fluctuations due 
to Alfven waves are damped; whereas in the region k > k c , Alfvenic fluctuations are faster than inertial 
waves and an equipartition of energy is still observed. 

5.3 Spectral anisotropy 

This section is devoted to the characterization in spectral space of the anisotropy of rotating MHD 
turbulence at high magnetic Reynolds number. The anisotropy of the flow is described using different 
statistical quantities. 

5.3.1 Shebalin angles 



The Shebalin angle 9q is often used in the study of MHD turbulence (Vorobcv et al. 2005 Shebalin 



et al. 



1983) and is defined, for any spectral quantity Q(k,t), by 



tan2 ^E fc feglQ(M)l 2 
E*W(M)l a 
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Figure 10: (Color online) Shebalin angles versus dimensionless time t* = tu /l . (a) 6 u (t*). (b) 6b(t*). 



where k^ — + fc^ is the horizontal component of the wave vector fc, whereas k z is the vertical 

component (aligned with the rotation axis). For an isotropic quantity, 9q s» 55°. A Shebalin angle of 90° 
is characteristic of a quantity independent of the vertical direction. This quantity and similar indicators 
are used in Favier et aL| ( 20io| for quasistatic MHD turbulence. 

The Shebalin angles of the kinetic and magnetic energy spectra Ek{U, 9) and Ej^(k, 9), 9 U and 9b, are 
gathered in figures |lO|a) and (b) respectively. Without rotation, both velocity and fluctuating magnetic 
fields are nearly isotropic with a Shebalin angle slightly greater than the isotropic value. This reflects a 
slight accumulation of kinetic and magnetic energies on equatorial modes such that k _L _B , due to the 



damping of nonlinear transfers along the imposed magnetic field ( Oughton et al. 1994 Alexakis et al 



2007). Note that, in the non-rotating case, the anisotropy is moderate since the interaction parameter 



is relatively small compared to previous studies (Bigot et al. 2008). However, as the rotation increases 



[i.e. as the Elsasser number decreases), both fields become strongly anisotropic since the two Shebalin 
angles increase. A Shebalin angle close to 90° is characteristic of a flow which tends to be invariant in 
the vertical direction, i.e. in which the vertical gradients d/dz decrease strongly, or, equivalently, the 
spectral energy concentrates in equatorial modes such that k _L fi. Unlike the monotonous evolution 
of 9b with A, after a first increase with A from oo to 0.5, 9 U decreases at very low Rossby numbers, as 
illustrated by the curves at A = 2 and A = 0.5 on figure 10 [a). This effect has already been observed 



in hydrodynamic rotating turbulence at very low Rossby number (Cambon et al. 1997 Morinishi et al 



2001 1. Very high rotation rate leads to an almost complete inhibition of nonlinear kinetic energy cascade, 



so that the anisotropy is reduced. 

The growth of 9b in time from the isotropic value is faster than the growth of 9 U for all finite values of 
A, indicating a faster structuring of the magnetic field towards a quasi-twodimcnsional state. A possible 
explanation lies in the nature of the phenomenon responsible for the growth of magnetic energy. Given 
that Alfven waves are linear, the anisotropic growth of the magnetic energy is much faster than the slow 
nonlinear anisotropic energy transfer due to the Coriolis force. Note that the two-dimensionalization 
of both the velocity and the magnetic fields is only partial since the Shebalin angles are far from the 
extremal value of 90°, characterizing a field invariant in the vertical direction. 



5.3.2 Angle dependent spectra of two-point correlations 

The Shebalin angle is a simple scalar characterizing the global anisotropy of the flow because of the 



integration over all wave vectors in equation ( 24 1 . In order to describe the anisotropy over all scales of the 



flow, it is necessary to introduce a scale- and angular-dependent quantity. The angular energy spectrum 
is a very useful quantity in the axisymmetric case. Contrary to classical energy spectra already presented 
in figures [9j the energy density spectrum e(k) is not spherically averaged, leading to E(k) — J e(fe)d 3 fc, 
but the angular dependence on the polar angle 9 is conserved : 



E K (k,9) 



8+A8/2 



cos9d9 



-A8/2 



E 



u t (k,9)u*(k,9) 



(25) 



k-Ak/2<\k\<k+Ak/2 
8-A8/2<6<8+Ae/2 
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Figure 11: (Color online) Angular energy spectra at t* w 5. Red color (or dotted lines) corresponds to 
the pole (i.e. modes such that 9 « 0) whereas blue color (or continuous lines) corresponds to the equator 
(i.e. modes such that 9 w tt/2). The lines without symbols correspond to the toroidal energy whereas the 
lines with symbols represent the poloidal energy, (a) Angular kinetic energy spectra without rotation, 
(b) Angular kinetic energy spectra with rotation at the same time, (c) Angular magnetic energy spectra 
without rotation, (d) Angular magnetic energy spectra with rotation at the same time. 



where Afc and A9 specify the discretization steps in Fourier space used for computing the anisotropic 
spectra (see figure [I] the blue region corresponds to the scales which contribute to E(k,9)). Here, 
Afc = 1 and A9 = n/10, these figures are necessarily linked to the DNS resolution, to ensure good 
statistical sampling. In the isotropic case, all angular spectra collapse. 

In addition to angular dependence, it is convenient to introduce the Craya-Herring frame (e^ , , e^ 3 ) ) 
(Sagaut and Cambon 2008) to geometrically decompose the velocity field. Because of the incompressibil- 
ity condition V-it = and because of V-6 = 0, the velocity and the fluctuating magnetic field spectral 
components along e^ 3 ) are identically zero since k-u = and k-b = 0. Any spectral component can 
therefore be decomposed as a toroidal contribution (along e^ 1 )) and a poloidal contribution (along e^). 
This representation is linked to the decomposition between shear (poloidal) and pseudo (toroidal) Alfven 



waves (Bigot et al. 20081 



The angular energy spectra are plotted in figure [IT] at t* « 5 with and without rotation, for A — > oo 
and A = 0.5. The equator (modes such that k _L fi) corresponds to blue solid lines whereas the pole 
(modes such that fe || O) corresponds to red dotted lines. For the sake of clarity, only polar and equatorial 
modes are represented (the spectra at intermediate angles < 9 < tt/2 lie monotonously between the 
spectra at 9 = and 9 — tt/2). The lines without symbols correspond to toroidal spectra whereas the 
lines with symbols correspond to poloidal spectra. In the following, we focus on the inertial and small 
scale dynamics (i.e. fc > 6) since the statistical sampling for computing the angular-dependent spectra 
at low wave numbers is too coarse. 

Let us first consider the velocity spectra. In absence of rotation, figure [TT| a) shows that the velocity 
field is roughly isotropic at all scales since all directional spectra almost collapse: there is equipartition 
of kinetic energy between polar and equatorial modes and equipartition of energy between poloidal and 
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toroidal components. One observes however a small increase of equatorial energy with respect to the 
polar one in the inertial range, which is consistent with a Shebalin angle slightly greater than its isotropic 
value, already seen on figure [lO{a) . 

With rotation, figureJTTJb) shows that the anisotropy of the velocity field is similar to the one observed 
in non-magnetized rotating flows. Without discriminating poloidal and toroidal contributions, one first 
observes a concentration of kinetic energy at the equator. Surprisingly, this anisotropy increases with 
the wave number, which is apparently inconsistent with a possible restoration of isotropy at small-scales, 
even if one accounts for the low Reynolds number of the simulations. By analogy with the Ozmidov scale 



in stratified fluids, some authors (Zeman 1994 Zhou 1995) introduce a cut-off wave number 



k Q = y — (26) 

where is the kinetic dissipation rate. When k > fcn, the effect of the Coriolis force may be neglected 
and a restoration of isotropy is expected. When k < fen, anisotropic effects are dominant. In figure [TT^ b), 
fen ~ 280 so that, based on this approach, all scales are dominated by rotation. Thus, higher Reynolds 
number and numerical resolution would be required to observe the restoration of isotropy at small scales. 

Let us now distinguish the toroidal kinetic energy spectrum (solid lines) from the poloidal kinetic 
energy spectrum (lines with symbols). This decomposition is of particular interest for the equatorial 
modes. At the pole, axisymmetry imposes the poloidal and toroidal components to be equal. One indeed 
observes equipartition between poloidal and toroidal energy for polar modes but not for equatorial modes, 
for which the toroidal energy is dominant at large scales whereas the poloidal energy is dominant at small 
scales. This kind of anisotropy has already been observed in MHD turbulence at very low magnetic 



Reynolds number using the quasi-static approximation (Favier et al. 2010). In that case, anisotropic 
Joule dissipation leads to a flow invariant in the direction of the imposed magnetic field. From that 
quasi- twodimensional state in which d/dz sa but u z is non necessarily zero, the non-linear cascade of 
energy can be very different depending on the component considered. The vertical component of velocity 
behaves like a passive scalar (with a classical cascade) whereas the horizontal component behaves like 
purely 2D turbulence and thus displays a weak direct cascade and a possible inverse cascade in the 
absence of forcing. This observation is consistent with the dominance of toroidal energy at small scales, 
with a horizontal toroidal mode, as shown on figure [T] However, in rotating turbulence, the mechanisms 
responsible for the quasi-twodimensionalization differ from linear dissipative effects of quasi-static MHD 
turbulence. The transition from 3D to quasi-2D is triggered by non-linear energy transfers and is thus 
less effective (Cambon 2001). The previous scenario for quasi-static MHD turbulence, valid in low Rm 



turbulence, is nonetheless very promising in the present case. 

Let us now discuss the magnetic energy spectra depicted on figures |TT|c) and (d) . Without rotation, 
the equipartition of energy is observed and the anisotropy of the velocity and magnetic fields are very 
similar. With rotation, one observes again the attenuation of the magnetic energy at large scales, and 
equipartition of energy at small scales. Thus, the anisotropy of the small-scale velocity field imposes the 
anisotropy of the small-scale magnetic field, with dominant equatorial energy. The fluctuating magnetic 
field is therefore nearly invariant in the vertical direction, which is also consistent with the previous 
observation concerning the increase of the Shebalin angle 0^. Note however that, focusing on equatorial 
modes (blue lines) , the fluctuating magnetic field is characterized by a dominant poloidal energy at large 
scales and a dominant toroidal energy at small scales, which is the opposite of the velocity field. 

In the rotating cases, both velocity and magnetic energies are concentrated inside equatorial modes. 
However, for a given wavevector k _L 17, the Fourier component can be vertical (i.e. all the energy is 
therefore poloidal) or horizontal (i.e. all the energy is toroidal, see figure [l]), or in between (with mixed 
poloidal and toroidal contributions). The previous analysis seems to indicate that the small scale velocity 
field is dominated by its vertical component (the poloidal energy being slightly dominant on figure [TT^ b) ) 
whereas the small scale magnetic fluctuations are dominated by their horizontal components (since the 
toroidal energy is dominant on figure [Tljd)) . This conclusion is confirmed when looking at the spherically 
averaged spectra, computed component by component. We separate the vertical energy spectrum 

Ei(k)= J2 Mwm (27) 

fc-Afc<|fc|<fc+Afc 

from the horizontal energy spectrum 

E\{k) = \ J2 [fii(fcK(fc)+fi2(fc)fi5(fc)] (28) 



k-Ak<\k\<k+Ak 
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Figure 12: (Color online) Scale-by-scale ratio between horizontal and vertical energies at time t* ~ 5. 
The lines represent the kinetic energy spectra whereas the lines with symbols represent magnetic energy 
spectra. 

where the index 1 and 2 refer to the directions perpendicular to Bq and £1, and the index 3 refers to 
the vertical direction. Similar definitions apply for the magnetic fluctuations. Figure [12] presents the 
scale-by-scale ratio between horizontal and vertical energy spectra of magnetic and velocity fluctuations. 
Without rotation (blue color), E\(k) / E v K (k) and E\j{k) / 'E v M (k) are about unity at all scales. For A = 0.5 
(red color), one indeed observes that the vertical kinetic energy is dominant in the inertial range (i.e. 
E\jE v K < 1) whereas the horizontal magnetic energy is dominant (i.e. E^/E^ > 1). This confirms 
quantitatively our previous observations based on the poloidal-toroidal decomposition. 

In summary, the anisotropy of the velocity field for rotating MHD turbulence at high magnetic 
Reynolds number is very similar to the anisotropy observed in non-magnetized rotating flows. The 
non-linear angular transfers due to rotation accumulate the energy toward the equatorial plane. This 
quasi-twodimensionalization leads to a departure from the equipartition between toroidal and poloidal 
energies. The fluctuating magnetic field is damped at large scales due to inertial waves, and is strongly 
anisotropic at small scales. The small scale velocity field is dominated by vertical motion, whereas the 
small scale magnetic field is mostly horizontal. 

6 Conclusion 

We present in this paper data from direct numerical simulations of homogeneous incompressible turbu- 
lence submitted to both Coriolis and Lorentz forces. This type of flows are of geophysical and astrophysical 
interests but are however characterized by numerous dimensionless parameters. We therefore focus on 
the large magnetic Reynolds number, moderate interaction parameter and small Rossby number regime. 

We have focused here on the study of the turbulent induction. In that case, the equipartition between 
kinetic and magnetic energy due to Alfven waves is broken by inertial waves introducing a separation 
scale, which depend on the Lehnert number. This departure from equipartition leads to a reduction of the 
large scales magnetic energy in aid of the kinetic energy. Concerning anisotropy, which is a key element 
in these flows, the velocity field is very similar to the well documented rotating hydrodynamic turbulence 
(because the anisotropy due to the imposed magnetic field is small) , with a concentration of energy in the 
modes perpendicular to the rotation axis, which corresponds to the two-dimensional manifold in physical 
space. The magnetic fluctuations display the same kind of angular anisotropy as the velocity field (i.e. 
dominant equatorial energy at small scales) , but with a different repartition between poloidal and toroidal 
components: the small-scale velocity field is dominated by poloidal modes (as in hydrodynamic rotating 
turbulence) whereas the small-scale magnetic field is dominated by toroidal modes. 

Many of the features of rotating MHD turbulence can be explained thanks to the linear analysis. 
Some basic properties of magneto-inertial waves (damping of the magnetic energy, misalignment between 
u and b) are of interest to study the fully nonlinear turbulent regime. 
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The small magnetic Reynolds is also of interest but the dynamics are in that case much more complex 
since both Coriolis force and Joule dissipation induce strong anisotropy on the flow. The angular transfer 
of energy due to Coriolis force compete with the anisotropic Joule dissipation, for example considering 
the quasi-static approximation. These two effects tend to concentrate kinetic energy on equatorial modes 
(i.e. modes such that k _L $7) and interesting dynamics could arise depending on the Elsasser number. 
The weak turbulence state (i.e. B ^> Uo), considered in previous studies (Bigot et al. 2008), could also 
be of interest. 

An other limitation of the present paper comes from the direction of the vector rotation O with respect 
to the imposed magnetic field Bo. We focus here on the axisymmetric case and the perpendicular case 
(being the most plausible configuration inside the Earth's core) is postponed to future studies. In that 
case, the competition between Alfven and inertial waves is much more complex, depending not only on 
the scale and on the polar angle 9 between k and Jl (being vertical) but also on the azimuthal angle 
between k and Bo (being horizontal). 

Another perspective concerns the forcing by physical instabilities. Decaying turbulence is chosen here 
not to alter the development of anisotropy. It is however possible to inject energy by mean of precessional 
instability (for instance by extending the theoretical study by Salhi and Cambon (2009) to DNS) or by 
unstable stratification using the homogeneous approach of Borue and Orszag (1997). 

Finally, initial conditions including initial magnetic fluctuations (with or without imposed magnetic 
field) is the next step to consider the effect of the Coriolis force on developed MHD turbulence. In that 
case, the simplified linear solutions (14) and (15) will depend on the initial condition for the magnetic 
mode, and a different dynamical regime is expected both from linear and nonlinear approaches. 

The authors would like to acknowledge support of the IDRIS (Institut du Developpement et des 
Ressources en Informatique Scientifique) for computational time on NEC SX-8 under the project 082206. 
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